Continuous weak measurement of the macroscopic quantum coherent oscillations 
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The problem of continuous quantum measurement of coherent oscillations in an individual quantum 
two-state system is studied for a generic model of the measuring device. It is shown that for 
a symmetric detector, the signal-to-noise ratio of the measurement, defined as the ratio of the 
amplitude of the oscillation line in the output spectrum to background noise, is independent of the 
coupling strength between oscillations and the detector, and is equal to (h/e) 2 , where e is the detector 
energy sensitivity. The fundamental quantum limit of 4 imposed by this result on the signal-to-noise 
ratio of the measurement with an "ideal" quantum-limited detector reflects the general tendency of 
a quantum measurement to localize the system in one of the eigenstates of the measured observable. 
These results are applied to specific measurements of the quantum oscillations of magnetic flux with 
a dc SQUID, and oscillations of charge measured with a Cooper-pair electrometer. They are also 
used to calculate the energy sensitivity of a quantum point contact as detector. 



I. INTRODUCTION 

Quantum superposition of macroscopically distinct 
states is one of the most characteristic features of quan- 
tum behavior at the macroscopic level. In the "meso- 
scopic" regime, when the states involved in the super- 
position correspond to a collective motion of a num- 
ber of particles that is larger than one but not quite 
macroscopic, the superposition of states has been demon, 
strated for photons in a high-quality microwave cavity 
and for the center-of-mass motion of large molecules 
Among the most basic dynamic manifestations of quan- 
tum superposition of states are quantum coherent oscil- 
lations between the two basis states of a two-state sys- 
tem. However, while the macroscopic quantum phenom- 
ena brought about by the incoherent quantum tunnel- 
ing are by now commonly found in a variety of systems 
ranging, from mesosccwpAc tunnel junctions in the engimes 
of fluxO'Q ancLchargeB'El dynamics, to molecularQ'Ei and 
nano-magnetsQ, the situation with experimental obser- 
vation of macroscopic quantum coherent (MQC) oscilla- 
tions remains much more uncertain. Claim of the obser- 
vation of the MQC oscillations in molecular magnetsEJ re- 
main highly ccptroversialrTri. Remarkable experimental 
demonstrationEll of the MQC oscillations in the charge- 
dynamics regime of a small Josephson junction is open 
to criticism that the two charge states of the observed 
quantum superposition differ in charge only by the the 
charge of one Cooper pair, and the oscillations between 
these two states can not be interpreted as macroscopic. 
Although this criticism is not fully justified, since the 
charge dynamics of a small Josephson junction is just 
another representation of its flux dynamics, which is-.the 
paradigm of the "macroscopic" quantum dynamicsU, it 
does not allow to consider the question of MQC oscilla- 
tions to be completely settled. 

Recently, macroscopic quantum dynamics has at- 
tracted renewed attention as the possible basis for de- 



velopment of scalable quantum logic circuits for quan- 
tum computation. In this context, macroscopic quan- 
tum two-state system plays the role of a qubit, an ele- 
mentary building block of a quantum computer. Several 
variants pf .qubits and quantum logic gates have been 
proposedtJlia that are based on the macroscopic quan- 
tum dynamics of Josephson junctions. Many character- 
istics of the macroscopic qubits compare favorably with 
those of the microscopic qubits: they are insensitive to 
disorder at the microscopic level and offer much larger 
freedom in design and fabrication of complex systems of 
qubits. The price of these advantages is the problem of 
the environment-induced decoherence, which is typically 
much more serious for the macroscopic than microscopic 
quantum systems. Suppression of decoherence to the ac- 
ceptable level requires thorough isolation of the qubit 
from its environment, the condition that typically limits 
the ability to control the qubit dynamics. 

This trade-off between the external control and de- 
coherence in a qubit has a fundamental aspect related 
to measurement. Even in the case of a perfect set-up, 
the measurement necessary to observe the MQC oscil- 
lations perturbs the system by projecting its state on 
the eigenstates of the measured observable, and there- 
fore presents an unavoidable source of decoherence. The 
intensity of this measurement-induced decoherence in- 
creases with increasing coupling strength between the 
detector and the oscillations, and more efficient measure- 
ment leads to stronger decoherence. First approaches to 
the problem of measurement ofJd^MQC oscillations in 
an individual two-states systemtlrEj~LiJ suggested for os- 
cillations of magnetic flux in SQUIDs, considered only 
the conventional limit of strong or "projective" quantum 
measurements, in which the detector-oscillation coupling 
is strong. In this case, the measurement leads to rapid lo- 
calization of the measured observable (flux) in one of its 
eigenstates, and suppresses the oscillation. This means 
that the time evolution of oscillations can be studied with 
strong measurements only if the detector can be switched 
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on and off on the time scale shorter that the oscillations 
period, and only in the "ensemble" of measurements, i.e. 
when the experiment is repeated many times with the 
same initial conditions. The measurement cycle consists 
then of preparation of the initial state of the system fol- 
lowed by its free evolution and the subsequent measure- 
ment. The information about dynamics of the oscillation 
is contained in the probability distribution of the mea- 
surement outcomes. Since the oscillation frequency is 
limited from below by several factors including the de- 
coherence rate and temperature, the need to switch the 
detector on and off rapidly in this approach presents at 
the very least a serious technical challenge. 

The goal of this work is to study quantitatively a new 



tions that 
, in which 



approach to measurement of the MQC op, 
is based on weak quantum measurementsca 
the dynamic interaction between the detector and the 
measured system is weak and does not establish perfect 
correlation between their states. Such a weak measure- 
ment provides only limited information about the system 
but, in contrast to strong measurements, perturbs the 
system only slightly and can be performed continuously. 
In this work, the process of continuous weak measure- 
ment of the MQC oscillations in an individual two-state 
system is considered quantitatively. Recent resultsEa for 
continuous measurements of electron oscillations in cou- 
pled quantum dots by a quantum point-contact are re- 
formulated within a generic detector model and applied 
explicitly to the MQC oscillations of flux measured by 
a dc SQUID and oscillations of charge measured with a 
Cooper-pair electrometer. It is shown that the signal-to- 
noise ratio of the measurement, defined as the ratio of the 
amplitude of the oscillation line in the output spectrum 
of the detector to the background noise, is fundamentally 
limited by the trade-off between the acquisition of infor- 
mation and dephasing due to detector backaction on the 
oscillations. This limitation is the least restrictive for a 
symmetric detector, for which the signal-to-noise ratio 
can be expressed as (U/e) 2 , where e is the detector en- 
ergy sensitivity. Since the energy sensitivity e is limited 
for regular (non-QND) quantum measurements by h/2, 
the signal-to-noise ratio of the continuous weak measure- 
ment of the quantum coherent oscillations is limited by 
4. This limit reflects the fundamental tendency of quan- 
tum measurement to localize the system in one of the 
eigenstates of the measured observable. As a spin-off, 
the established relation between the signal-to-noise ratio 
and energy sensitivity is used to demonstrate that the 
quantum point contact is the quantum-limited detector 
with energy sensitivity e = h/2. 



II. CONTINUOS MEASUREMENT OF THE MQC 
OSCILLATIONS WITH A LINEAR DETECTOR 



We consider the MQC oscillations in an individual two- 
states system, with the two basis states separated in en- 



ergy by e and coupled by the tunneling amplitude — A/ 2. 
The basis states are chosen to coincide with the eigen- 
states of an oscillating variable x, for instance, magnetic 
flux in a SQUID loop. In this basis, x — xq(t z /2, where 
xq is the difference between the values of x in the two 
states of the system, and a z is the Pauli matrix. In the 
simplest measurement scheme considered in this work, 
the detector measures directly the oscillating variable x, 
and therefore is coupled to x. This means that the Hamil- 
tonian describing the measurement set-up (Fig. 1) is: 



H=- l -{ea z 



Aa x + cr z f) + H 



(1) 



where Hq is the Hamiltonian of the detector, and / is the 
detector operator that couples it to the oscillations. For 
convenience, the amplitude xq of the oscillations and the 
coupling strength are included in /. 



- A/2 
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FIG. 1. The diagram of the MQC oscillations in an individ- 
ual two-state system measured by a linear detector. Detec- 
tor with the Hamiltonian Ho is weakly coupled through the 
operator / to the oscillating coordinate of the system. The 
detector output o(t) consists of noise and linear response to 
the input oscillations. 

We assume that the characteristic response time of the 
detector is short, e.g., much shorter than the period of 
measured oscillations, and that the detector operates in 
the linear regime. These two assumptions mean that the 
detector output o(t) can be written as 



o(t) = q(t) + -a z (t) 



(2) 



where q(t) is the noise part of the output, and A is a re- 
sponse coefficient. The linearity assumption implies that 
the coupling —<r z f/2 of the detector to the oscillations 
is sufficiently weak so that the detector response can be 
described in the linear-response approximation: 



A = 



dre 1 ' 



[Q(r),f])o- 



(3) 



The average (. . .) in eq. (||) is taken over the stationary 
density matrix of the detector. Since the response coef- 
ficient A is independent of frequency oj, eq. (||) can be 
written in terms of the q—f correlator as 



(q(t + T)f(t)) = 2*S q/ 8(T-0), S qf 



iX 



■in 



, (4) 



where the infinitesimal shift in the argument of the 5- 
function represents small but finite response time of the 
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detector, and is needed to resolve the ambiguity in eq. 
(||) . Parameter a in eq. (Q) is introduced to represent the 
real part of the q—f correlator that is not determined by 
the response coefficient A. 

The role of detector in a quantum measurement is to 
convert the quantum input signal, in our case, the oscil- 
lations x(t) = xocr z (t)/2, into the output signal that is 
already classical and can be dealt with (e.g., monitored or 
recordered) without "fundamental" problems. Condition 
of the classical behavior of the detector output requires 
the output spectral density to be much larger than the 
spectral density of the zero-point fluctuations in the rel- 
evant range of low frequencies of the input signal. For 
a detector with a short internal time scale this means 
that the noise q(t) is ^-correlated on the time scale of the 
input signal: 



(q(t + T)q(t)) =2irS 9 8(T). 



(5) 



Here S q is the constant low-frequency part of the spectral 
density S q (ui) of the detector output noise. 

For a quantum-limited detector at small temperature 
T — > 0, eqs. (Q) and impose a constraint on the spec- 
tral density Sf of the coupling operator /. Indeed, eqs. 
(f|) and (g) can be written explicitly in the basis of the 
energy eigenstates \k) of the detector and give the follow- 
ing expression for the spectral density S q : 



d£ k de k >v(e k )v(e k >)p k {k\q\k'){k'\q\k)5{£ k 



(6) 



Expression for the correlation amplitude S q f is similar, 
with the matrix element (k'\q\k) replaced by the matrix 
element of /. In these expressions, e k is the energy of the 
state \k), p k is the probability to be in this state, and v is 
the state energy density. Since S q is independent of the 
frequency u>, eq. (JsJ) is satisfied when the matrix elements 
of q and the density of states are constant, and eq. @ 
can then be written as 



S„ 



2,, 2 



It should be noted that the constant matrix elements 
(q) and (/) are off-diagonal in the fc-basis and can be 
imaginary. Following the same steps for Sf we express 
this spectral density in terms of S q and S q f. 



S f = \(f)\ 2 v 2 = \S qf \ 2 /S q 



(7) 



Equations (fjj) and (Q) relate the backaction noise of the 
detector determined by the spectral density Sf to its re- 
sponse coefficient and the output noise. Conceptually, 
such a relation resembles the fluctuation-dissipation the- 
orem that links response of the system to its equilib- 
rium fluctuations, but it does not have the status of a 
"theorem" . It is obvious from the derivation above that 
eq. (0) is not necessarily valid for an arbitrary system 
playing the role of detector in a quantum measurement. 



Nevertheless, it holds for several of the "standard" de- 
tectors: quantum point contact, resistively-shunted dc 
SQUID, and Cooper-pair electrometer, considered later 
in this work. 

Making use of q- and /- correlators, we can calculate 
the spectral density of the detector output o(t) in the pro- 
cess of continuos measurement of the MQC oscillations. 
From eq. (||), the correlation function of o(t) is: 

A 2 

K (t) = 2irS q S(T) + —Ty{p<j z <j z (t)} , (8) 

where p is the stationary density matrix of the two-state 
system established as a result of the interaction with the 
detector. Averaging the Heisenberg equation of motion 
of the operator <j z (t) over the 5-correlated backaction 
noise / of the detector, we get the set of equations for 
the time evolution of the matrix elements try of c z (r): 



<Tii = AImcri2 , <Ti2 = (is — r)<ri2 — iA(7n 
and (T22 = — <Tn, with the rate 



/ 



(9) 



(10) 



describing the backaction dephasing of the oscillations 
by the detector. The density matrix p of the two-state 
system satisfies the same set of equations (||), except for 
the normalization, pn + P22 — U and its stationary value 
is p = 1/2. Solving eqs. (||) with the initial condition 
0z(O) = o z and averaging a z a z (r) over p — 1/2 we find 
the spectral density S a (u) = (l/2n) K (r)e iUT dr. 
Under the conditions of "resonance", e = 0, when the 
oscillation amplitude is maximum, we get: 



S (uj) = S q + 



FX 2 



Air (w 2 -A 2 ) 2 



(11) 



When e 7^ 0, it is convenient to calculate the spectrum 
numerically from eq. (||). The spectrum in this case is 
plotted in Fig. 2 for several values of e and the dephas- 
ing rate T. For weak dephasing, r < A, the spectrum 
consists of a zero-frequency Lorentzian that vanishes at 
s = and grows with increasing |e|, and a peak at the 
oscillation frequency fl = (A 2 + E 2 ) 1 / 2 . The peak at zero 
frequency reflects the incoherent transitions with a small 
rate of order T between the states of the two-state sys- 
tem. The high-frequency peak of the MQC oscillations 
also has the width T. While this width can be small for 
sufficiently weak dot-contact coupling, the height of the 
oscillation peak cannot be arbitrarily large in comparison 



to the background noise spectral density S q . At 



0, 



when the amplitude of the oscillations is maximum, the 
peak height is S max — A 2 /47rr. Even in this case, the 
ratio of the peak height to the background is limited: 



A 2 



4A 2 



S„ 



ATt 2 S f S, 



A 2 + a 2 



< 4. 



(12) 



This limitation is universal, e.g., independent of the 
coupling strength between the detector and oscillations, 
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and reflects quantitatively the interplay between mea- 
surement of the MQC oscillations and their backaction 
dephasing. The fact that the height of the spectral line 
of the oscillations can not be much larger than the noise 
background means that, in the time domain, the oscilla- 
tions are drowned in the shot noise. 




co/A 
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FIG. 2. Spectral density S (u) of the output of the sym- 
metric detector with minimal backaction dephasing (a — 0) 
measuring the MQC oscillations for several values of (a) 
the energy bias e reflected in the oscillation frequency 
Q = (A 2 + e 2 ) 1/2 , and (b) the rate T @ of the measure- 
ment-induced dephasing. 



When the backaction dephasing rate T increases, the 
oscillation line broadens towards the lower frequencies, 
and eventually turns into the growing spectral peak at 
zero frequency associated with the incoherent jumps be- 
tween the two basis states of the two-state system. At 
large T, when the coherent oscillations are suppressed, 
the rate of incoherent tunneling decreases with increas- 
ing r. For instance, at T 3> fi, the tunneling rate is 
7 = A 2 /2r , and the spectral density of the detec- 
tor response at low frequencies to ~ 7 has the standard 



Lorentzian form, S (ui) — S q = 2j\ 2 /47r(47 2 + c 



Sup- 



pression of the tunneling rate 7 with increasing dephas- 
ing rate T is an example of the generic "Quantum Zeno 
Effect" in which quantum measurement suppresses the 



decay rate of a metastable state. In the context of search 
for the macroscopic quantum coherent oscillations, the 
Lorentzian spectral density has been observed and used 
for measuring the tunnelipff rate of incoherent quantum 
flux tunneling in SQUIDaSa. 

The maximum signal-to- noise ratio S max / S q (|lj) is at- 
tained if the fundamental backaction of the detector is the 
only mechanism of dephasing of the coherent oscillations. 
We now discuss briefly the effect of a weak additional de- 
phasing and energy-relaxation on the spectral density of 
the oscillations. The efect of such a weak relaxation is no- 
ticeable if the backaction dephasing is also weak, r <C A. 
Energy relaxation arises typically due to interaction with 
some external system ( "reservoir" ) that is in equilibrium 
at temperature T. The interaction term in the Hamilto- 
nian can be written similarly to the interaction with the 
detector ([[]) as 



H c — —(T z fr 



(13) 



where f r is the reservoir force acting on the system. Un- 
der the assumption of the frequency-independent relax- 
ation rate, the standard free equilibrium correlator of this 
force is (see, e.g., Ell): 



(fr(t)f r (t + r))=a 



we" 



7T 1 - e""/ T 



(14) 



where parameter a characterizes the relaxation strength. 
Comparison of this correlator with the (^-correlated back- 
action noise of the detector shows that the detector (|l|) 
is acting effectively as a reservoir with temperature that 
is much larger that the energies of the two-state system. 

Energy relaxation makes the stationary average value 
of a z non- vanishing, and the output correlation function 
should now be calculated as 

\2 

K (t) = K q (r) + — [(<r t a t (r) + <j z (t)<j z ) - 2(a z ) 2 } . 

o 

(15) 

For weak coupling, it is convenient to find the time evo- 
lution of ct z (t) in the basis of eigenstates of the two-state 
Hamiltonian. In this basis, the Hamiltonian including 
interaction with the reservoir (and omitting temporarily 
the detector) is : 



(16) 



H r is the Hamiltonian of the reservoir. Heisenberg equa- 
tions of motion that follow from the Hamiltonian ( |l6| ) 
are: 



fr = i[H r ,f r ] , 



ea z )[f r ,H r 



& z = 2f r Aa y , &± = + 2sf r /n)a± =F iAa z /fl , 
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where a± = (a x ± ia y )/2. Integration of the first two of 
these equations gives in the first order in coupling (|l~3| ) 
to the two-state system: 

f r (t) = f<?)(t) - iA fdr |V[/(°)(t),/(°>(t')K(t') , 

where fl°^ is the free part of the fluctuating reservoir 
force in absence of coupling. Solving the second pair 
of the Heisenberg equations up to the second order in 
coupling, making the rotating-wave approximation, and 
tracing out the reservoir degrees of freedom with the help 
of the correlator (|i~4|), we get a set of equations for the 
evolution of the matrix elements Sij of the operator of 
the oscillating variable (given by & z (t) in the original 
"position" basis) in the eigenstate basis: 



%(r) = r e [ 



si2(r) = (is - r )si 



S22J 



(17) 



Initial conditions for these equations are: sn = — S22 — 
e/f2, and S12 = — A/O. The characteristic energy- 
relaxation rate in eq. (O) is T e = 2aA 2 /f2, and the total 
dephasing rate is 



2 n 



In eqs. ( \L7\) , we also added the detector dephasing terms 
from cq. (ft) "rotated" from the position basis into the 
eigenstate basis. 

The density matrix r of the two-state system in the 
basis of eigenstates satisfies similar equations, and the 
stationary values of its matrix elements are ri2 = and 
r n = (r t + r e )/2r t , where 

r t ee r e coth(f7/2T) + rA 2 /n 2 . 

Using these relations, the definition (|l^), and the evolu- 
tion equations ( |l7| ) we find the spectral density: 



S (lo) = S q + 



X 2 



2e 2 T t 



E 



(w ± 



v 2 
1 



(18) 



As in the case without energy relaxation, the spectral 
density consists of a zero-frequency Lorentzian of width 
T t and peaks at ±0 of width To due to coherent oscil- 
lations. For weak relaxation, the incoherent slow transi- 
tions giving rise to the low-frequency noise are the tran- 
sitions between the two energy eigenstates of the system. 
The height of the oscillation peak is suppressed in pres- 
ence of additional energy relaxation that contributes to 
the dephasing rate T e , and the relative magnitude of the 
peak, Smax/Sq is smaller than its value without the re- 
laxation. 



III. RELATION TO ENERGY SENSITIVITY 

The detector characteristics for measurement of the 
quantum coherent oscillations in a two-state system con- 
sidered in the previous Section are related to another de- 
tector characteristic that is used for meapir|eHients of the 
harmonic signals - see, e.g., examples inE3~Eil, and some- 
times is loosely referred to as "energy sensitivity". It is 
defined by considering the detector measuring a harmonic 
oscillator with a frequency u>q and a small relaxation rate 
7 < wo. In this Section, we establish the quantitative re- 
lation between the signal-to-noise ratio Smax I Sq for mea- 
surements of the two-state systems discussed above with 
the energy sensitivity used in the literature for measure- 
ments of harmonic signals. 

The Hamiltonian of the damped harmonic oscillator 
attached to a detector is obtained by replacing the two- 
state part of eqs. ([!]) and ( [l3|) with the corresponding 
oscillator terms: 



H= y(i 2 +c V) 



x(f + f r )+H + H r 



(19) 



where M is the mass of the oscillator and x is the os- 
cillating coordinate. Due to linearity of the system, the 
Heisenberg e qua tion of motion for x that follow from the 
Hamiltonian ( J19| ) (see, e.g. j£3) coincides with the classical 
equation of motion of the damped oscillator: 

x- 1 x + co 2 x = -L(f r (t) + f(t)). (20) 

As in the previous Section, the random forces f r (t) and 
f(t) are produced, respectively, by the reservoir respon- 
sible for the energy relaxation of the oscillator and by 
the detector. (Although the operators f r and / in eqs. 
(19) and (pp|), as well as the transfer coefficient A in eq. 
(21) below, differ from the corresponding quantities used 
in Section 2 by a normalization factor, this distinction is 
not made explicit. This should not lead to any confusion, 
since the normalization factor drops out of all expression 
which are compared between the two Sections.) In gen- 
eral, the right-hand-side of eq. (^0|) should also contain an 
external perturbation that creates a "signal" component 
of the oscillations x(t). However, for the discussion of the 
detector sensitivity, it is appropriate to treat the equilib- 
rium fluctuations of the oscillator driven by the reservoir 
noise f r (t) (e.g., the zero-point oscillations at vanishing 
temperature T) as part of the signal. This allows us not 
to include additional signal terms in eq. (p0|). 

The sensitivity of the detector is characterized by the 
detector noise contribution to the spectral density of the 
output S reduced to the detector input. Similarly to 
eq. (^|) the output of the detector measuring harmonic 
oscillator is: 



y(t) = q(t) + Xx(t) 



(21) 



This equation implies that the detector contribution to 
the output noise comes from the two sources: direct out- 
put noise q(t), and effect of the backaction noise f(t) on 
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the oscillator coordinate x{t). Introducing the dynamic 
response function G(t) of the oscillator: 



and 



X{t)= drG(T)(f r (t-T) + f(t-T)), (22) 

Jo 

with G(t) = for t < 0, we see from eq. (|2(]) that 

G{uo) = / dre-^G(r) = * ^— . (23) 



■ IJLO 



In terms of the response function, x(co) = G(u>)(f r (uj) + 

/M)- 

Since the detector noise /(<) is uncorrelated with the 
reservoir force f r (t) , or in general, any other signal com- 
ponent of x(t), we see from eqs. (pi]), (p2]), and ( p3| ) that 
the spectral density of the detector output consists of the 
two additive components: signal, i.e., equilibrium spec- 
tral density of the harmonic oscillations transformed to 
the output, and the detector noise Sn(lo): 

S (u) = X 2 \G{Lo)\ 2 S r (Lo) + S n {lo) , 



SVH =S q + X 2 \G(u))\ 2 S f + 2\ReG(uj)Sq f . (24) 

Here S q f is the symmetrized correlator of the detector 
output and input noises, S q f — (S q f + Sf q )/2 — a/4ir 
(with a defined in eq. (^)), and S r (to) is the spectral 
density of the reservoir force f r : 



S r (co) = ( 7 wM/27r)coth(?iw/2T). 



(25) 



The noise properties of the detector are better charac- 
terized if its contribution to the output noise is reduced 
to the input, i.e. instead of Sn(u) (^4|) we consider the 
quantity F(ui) = Sn(co)/ X 2 \G(lo)\ 2 . For weak damping, 
7 <C loo, when 



|GM| 2 



1 



4 W 2 M 2 ((^~^ ) 2 + 7 2 /4) ' 
we get 

F(lo) = S f + S q {2uo M/\) 2 {{uj - loo) 2 + 7 2 /4)- 

S q f (4woM/A) (lo - uo) ■ 

The three terms in this equation scale differently with the 
strength of the detector-oscillator coupling, since A and 
S q f are proportional to the first power, while Sf is pro- 
portional to the second power of the coupling strength. 
F(lo) can be minimized with respect to the coupling 
strength and also with respect to the small detuning 
lo — loq between the signal frequency and the oscillator 
frequency. The minimum is reached when 



LO 



XS, 



if 



-fco MS q 



Sq 



So — (SqSf 



dqf) 



and is equal to 



Fmin = 2jLO M 



Sq 

T 



(26) 



It is convenient to normalize the reduced noise F in 
such a way that it can be directly compared to the equi- 
librium fluctuations in the oscillator driven by S r - Since 
at this stage we can already neglect small difference be- 
tween the oscillator frequency and signal frequency, both 
the minimum noise f m j„ (Eq) and the equilibrium spec- 
tral density S r (p5j ) are proportional to the oscillator pa- 
rameters jlooM. This means that with appropriate nor- 
malization we can define F m i n directly in terms of the 
number of quanta added to the signal at frequency lo by 
the detector. This is achieved by introducing the energy 
sensitivity as 

Equations (Q) and (Q) of the previous Section show that 
for a quantum-limited detector 



e = h/2. 



(28) 



(Note that in this Chapter, Ti is shown only in some of the 
final results.) Equation (|28| ) agrees with the conclusion 
of a general theory of quantum linear amplifiers, accord- 
ing to which a phase-insensitive linear amplifier adds at 
least half-a-quantum of noise to the amplified signal - 
see, e.g.,c3. This result is related to eq. (ref39) since de- 
tector in the measurement process plays the role of an 
amplifier transforming weak quantum input signal into 
the classical output. 

Comparing eq. ( p7| ) for the energy sensitivity with eq. 
( |l2| ) for the signal-to-noise ratio of the measurement of 
the quantum coherent oscillations we see finally that 
these two quantities are closely related. When the out- 
put and input noises of the detector are uncorrelated, 
S q f = 0, the relation is simple: 



S a 



(h/e) 2 



(29) 



2tO MSq 



As will be clear from the examples of specific detector 
considered below, the situation with S q f = can be rea- 
sonably referred to as "symmetric detector" . As follows 
from eq. (ref40), the largest signal-to- noise ratio of 4 is 
obtained for such a symmetric detector in the quantum- 
limited regime with e — h/2. 

When the input-output correlation is non- vanishing, 
the signal-to-noise ratio is smaller that the value given 
by eq. (|2^), while the energy sensitivity e for measure- 
ment of the harmonic signal can still be made equal to 



G 



h/2 by optimizing the detuning between the signal and 
the oscillator. Another difference between the measure- 
ment of harmonic oscillator defining e and measurement 
of the two-state system, is that the minimum noise (^g|) in 
the oscillator measurement is reached only for optimum 
detector-oscillator coupling, while the maximum signal- 
to-noise ratio ([12]) for the quantum oscillation measure- 
ment is independent of the coupling strength to the de- 
tector as long as the coupling is weak. 



IV. ENERGY SENSITIVITY OF A QUANTUM 
POINT CONTACT 



One of the applications of the results obtained in the 
previous Section is the demonstration that a quantum 



y used as the detector of 
can reach the quantum- 



point contact that is fre> 
electric charge or voltagi 
limited regime with ultimate energy sensitivity (Gq) . The 
mechanism of operation of a quantum point contact as a 
detector utilizes modifications of the electron transmis- 
sion properties of the contact by the measured voltageE3. 
When the contact is biased with a large voltage V, 
changes in the electron transmission probability lead to 
changes in the current / flowing through the contact 
which serve as the measurement output. Fluctuations of 
the electric potential in the contact region due to the cur- 
rent flow produce the backaction dejphasing of the mea- 
sured object by the point contactE§E3. This dephasing 
was calculated-. fpe symmetric contacts within different 
approaches irtZMid. It is knownOLJ that in the case of 
measurement of a two-state system, when the coupling 
to the system is symmetric, quantum point contact is 
an ideal detector of the quantum coherent oscillations. 
Such a detector causes the minimum dephasing of the 
oscillations that is consistent with the information ac- 
quisition by the measurement. For asymmetric coupling, 
the dephasing by thejpaint contact is larger that the fun- 
damental minimuirou. 

To calculate the energy sensitivity of the point contact 
detector, we start with the standard Hamiltonian of a 
single- mode point contact. Including a weak additional 
scattering potential U(x) for the point contact electrons 
which is the input signal of the measurement we can write 
the Hamiltonian as 



H = £ kal k a ik 



ik 



ij kp 



JP • 



(30) 



The operators in this Hamiltonian represent point- 
contact electrons in the two scattering states i = 1, 2 (in- 
cident from the two contact electrodes) with momentum 
k, and Ujj = f dx ip% (x)U (x)tpj (x) are the matrix ele- 
ments of the potential U(x) in the basis of the scattering 
states. Here ipi( x ) is the wavefunction of the scattering 
state, and x is the coordinate along the point contact. 
Several assumptions are made about the contact. The 



bias energy eV is assumed to be much larger than tem- 
perature T, but much smaller than both the Fermi energy 
in the point contact and the inverse traversal time of the 
contact. This allows us to linearize the energy spectrum 
of the point-contact electrons: e k — v F k, where vp is 
the Fermi velocity, and neglect the momentum depen- 
dence of the matrix elements [Zy. The potential U[x) is 
also assumed to be sufficiently weak and can be treated 
as perturbation. In this regime, the point contact oper- 
ates as a linear detector, and the current response to the 
perturbation U can be calculated in the linear-response 
approximation. The last assumption is that the frequen- 
cies of the input signal are much smaller than eV, the 
fact that allows to treat U as the static perturbation. 

At frequencies much lower that both eV and inverse 
traversal time of the contact, the current is constant 
throughout the contact and the contact response can be 
calculated at any point x. We choose the origin of the 
coordinate x in such a way that the unperturbed scatter- 
ing potential is effectively symmetric, i.e., the reflection 
amplitudes for both scattering states are the same, and 
then take x to lie in the asymptotic region of the scatter- 
ing states. In this case, the standard expression for the 
current in terms of the electron operators ^{x): 

1=^ — * t T i 5-*), * (x = > Vifc (x )a ik , 

2m ox ox 

ik 



jives for the current operator at x: 

I = "V \D(a]i.a-\„ — a\ k a2 P ) 



evF 



kp 



(31) 



Here D and R = 1 — D are the transmission and reflec- 
tion probabilities of the point contact, L is a normaliza- 
tion length, and the variation of the momentum near the 
Fermi points (i.e., the difference between k and p) was 
neglected everywhere besides the phase factor in the sec- 
ond term. The reason for keeping this factor will become 
clear later. 

In the linear-response regime, the current response of 
the point contact is driven by the part of the perturbation 
U causing transitions between the two scattering states 
ipijt- As shown in the Appendix, the real part of the 
transition matrix element U12 is related to the change 
SD of the transmission probability of the contact: 



U 12 = 



v F 



5D 



L 2(DR) 1 / 2 



U21 = U 



12 ■ 



(32) 



The imaginary part of C/12, expressed through a dimen- 
sionless parameter u in eq. (Q) , does not affect the current 
/. Qualitatively, it characterizes the degree of asymmetry 
in the coupling of the quantum dots to the point contact; 
u = if the perturbation potential U (x) is applied sym- 
metrically with respect to the main scattering potential 
of the point contact. 
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In the measurement process, the perturbation U rep- 
resents the coupling operator between the contact and 
the measured system, whereas the current I is the mea- 
surement output. As follows from the discussion in the 
previous Section, the energy sensitivity of the contact is 
determined by the correlators of U and /, and by the 
response coefficient A. Using eqs. (|30|), (|3l|), and ( |32| ) we 
can evaluate the correlators directly. In the limit of large 
voltages, eV T, when the contact noise properties are 
dominated by the shot noise, we get: 



eV {8Df + u 2 



DR 

e 3 VDR 



(U(t)U(t + r)) = — 

(I(t + T)I(t)) 

71 

(U(t)I(t + r)) = ^(i5D + u) S(t - rj) . 

Z7T 



S(r) 
6(r) 



(33) 



The time delay rj = in the last of eqs. (|33f ) comes 

from the phase factor e~ l V°~ p " x \ kept in eq. (pi|), and is 
infinitesimally small for small traversal time of the con- 
tact. It is nevertheless important for correct calculation 
of the contact response A. From the U—I correlator and 
the standard expression for the linear response (^) we 
confirm that A is equal to the change of current through 
the contact due to change SD of the transmission coeffi- 
cient, A = e 2 VSD/n. 

The correlators ( |33| ) satisfy the general relations (0) 
and (^) , and therefore the energy sensitivity of the quan- 
tum point contact reaches the fundamental quantum 
limit (28). It should be noted, however, that this con- 
clusion is strictly valid only in the large-voltage limit 
eV ^> T. At finite temperature T, scattering of the 
point-contact electrons within the same direction of prop- 
agation (described by the terms Uu and U22 of the per- 
turbation U) creates additional contribution to the back- 
action noise and degrades the energy sensitivity. The 
magnitude of this effect depends on the magnitude of the 
"forward" scattering matrix elements U\\ and TJ22 rela- 
tive to the backscattcring matrix element U12 , increasing 
with intensity of the forward scattering but decreasing 
with T/eV ratio. 



V. FLUX AND CHARGE MQC OSCILLATIONS 

This Section provides specific examples of measure- 
ments of the macroscopic quantum coherent oscillations 
of magnetic flux and electric charge. It is shown that the 
typical detectors for the flux and charge measurements, 
dc SQUID and Cooper-pair electrometer, satisfy the gen- 
eral equations of Sections 2 and 3, and should be capable 
of reaching the fundamental limit of the signal-to-noise 
ratio for the continuous weak measurement of the MQC 
oscillations. 



A. Flux oscillations measured with a dc SQUID 

Typical set-up of a measurement of the MQC oscilla- 
tions of flux with a dc SQUID consists of a two-state flux 
system (rf SQUID with half of a magnetic flux quantum 
<&o = 7r?i/e induced in it by an external magnetic field) 
coupled inductively to a dc SQUID biased with an exter- 
nal current Iq and shunted by a resistor R (Fig. 3) . When 
the inductance of dc SQUID loop is small, the difference 
between the two Josephson phases ipi^ across the two 
junctions of the SQUID is directly linked to the flux $ 
induced in the dc SQUID by the flux oscillations: 

<4>\ - <P2 = 2vt$/<I>o = 6 . 

In this case, the SQUID is equivalent to a single Joseph- 
son junction, with the supercurrent in this junction mod- 
ulated by the flux $. The total amplitude of Cooper-pair 
tunneling in the SQUID is equal to a sum of the two in- 
dividual amplitudes of tunneling in the two SQUID junc- 
tion, and can be written as 

Ej/2 = 7 (+) (e)/4e, = he lB/2 + I 2 e- lB/2 , 

(34) 

where I1.2 are the critical currents of the two junctions. 
Coherent sum of the two tunneling amplitudes in ( |3~i| ) 
leads to modulation of the total supercurrent of the 
SQUID. 




FIG. 3. Schematic diagram of a continuous measurement 
of the MQC oscillations of magnetic flux with a dc SQUID. 
The two-state flux system is implemented as an rf SQUID 
with externally induced flux $o/2 and is coupled to the dc 
SQUID biased with a current Jo- The coupling strength is 
characterized by the variation <5$ of the flux through the 
dc SQUID due to oscillations. The voltage V across the dc 
SQUID is the measurement output. 

We consider the simplest regime of the dc SQUID dy- 
namics when the bias current Iq and associated average 
voltage Vo — Rio across the dc SQUID are sufficiently 
large, and the dc component of the Josephson current 
through it is small in comparison to Iq. In this regime, 
the Cooper-pair tunneling through the SQUID is ade- 
quately described by perturbation theory in the tunnel- 
ing amplitude (|3) and can be qualitatively interpreted 
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as incoherent jumps of individual Cooper pairs. The re- 
sistor R provides the dissipation mechanism that trans- 
forms reversible dissipationless Cooper-pair oscillations 
between the two electrodes of the SQUID into incoher- 
ent tunneling. 

Using the known results for the incoherent Cooper-pair 
tunnclingp3, we can calculate the rate of this tunneling 
and find all the detector characteristics of the dc SQUID. 
The SQUID is coupled to the oscillations by the operator 
of the current I- circulating in the SQUID loop multi- 
plied by the change 5§ of the flux in the loop induced by 
the flux oscillations. Changes in the flux through the dc 
SQUID change the total current 1+ through both SQUID 
junctions and create deviations V of the voltage across 
the SQUID from Vq, V — —RI+, that serve as the mea- 
surement output. As follows from the general discussion 
in Sec. 2, SQUID parameters important for measurement 
are the coefficient A of the transformation of the oscillat- 
ing flux into the voltage V, the spectral density Si of 
the circulating current I- that is responsible for back- 
action dephasing by the SQUID, the spectral density of 
the output voltage Sv, and the correlator Svi between 
V and These parameters can be found quantita- 
tively starting from the tunneling part Ht of the SQUID 
Hamiltonian that can be written as 



Hi 



E J A(2eV t+ V (t)) 



h.c. 



(35) 



with ip(t) being the random Josephson phase across the 
dc SQUID accumulated due to equilibrium voltage fluc- 
tuations produced by the resistor R. It is characterized 
by the correlator 



(¥>(*)¥>) 



u> 1 — e W / J 



(36) 



where p = R/Rq is the resistance R in units of the quan- 
tum resistance Rq — nh/Ae 2 , the average (. . .) is taken 
over equilibrium density matrix of the resistor R, and 
g(u>) describes the cut-off of the dissipation provided by 
R at some large frequency uj c associated with either fi- 
nite inductance of the SQUID or finite capacitance of its 
junctions, while = 1 at w < w c . 

The operators of the two currents I± that determine 
the SQUID parameters are: 



: [/ (±) (e)e l(2ey ° t+<p(t)) 



h.c 



(37) 



ee {he lQ ' 2 -I 2 e- lQ / 2 )/2. 



In the regime of the incoherent Cooper-pair tunneling 
the average dc current (/ + ) can be found treating the 
tunneling Ht as perturbation: 

(/+> = -*/ dt([I + ,H T (t)])=7r\I^(Q)\ 2 r/e, (38) 
Jo 

where 



T-R-e 

47T 



dte 



ileVat I 



[e tlp{t \e- lv ]) . 



(39) 



For example, for vanishing temperature T, and small bias 
voltages, 2eVo <C w c , the time r defined in (EEj) can be 
found from eqs. (|39|) and (36) to be (see, e.g.,H): 



(l/4^ c r(p))(2eUoM 



(40) 



When the resistance R is small, R <C Rq, t becomes 
independent of u) c , t = eR/2irV. In this case, all 
the SQUID characteristics, including the average current 
(|3^) can be obtained by direct time averaging of the 
classical Josephson oscillations in the SQUID. 

The noise spectral densities of the two currents, I±, 
are obtained by directly taking the average over equilib- 
rium density matrix of the resistor R. They vary with 
frequency on the scale of the Josephson frequency 2eVo, 
and arc constant at lu <C 2eVo- In this frequency range, 



dte Jielw ([e— ',e 



(41) 



Fluctuation-dissipation theorem relates the anticommu- 
tator [•••]+ in this equation to the commutator in eq. 
(E59|) and gives: 



5/ = |/(-)(e)|V, s v = i? 2 |/ (+) (e) 



V. 



(42) 



where r' = t coth(eVo/T). The correlation function Syi 
is found similarly: 



Svi =R[i (+ \e)]*i ( -\ey . 



(43) 



Comparison of the spectral density Sv ( ff2| ) and the 
average current ( |38| ) shows that Si + = Sv/R 2 = 
(e(I + )/ir) coth(eVo/2 n ), i.e. the noise of the current 1+ 
can indeed be interpreted as resulting from uncorrelated 
transitions of individual Cooper pairs. In particular, at 
T <C eVo, the noise is the shot noise of Cooper pairs. 

Finally, eq. ( |38|) gives the response coefficient of the 
SQUID 



A ee dV/d® = 2-KR(d\I (+) {Q)\ 2 /dQ)r. 



(44) 



(Note that as in eq. (|4^) for the backaction noise Si and 
also in eq. (^) for the correlator Syi, the factor (5$ is 
omitted from the definition of A.) 

For temperatures negligible on the scale of eVb, r' 
in the noise spectral densities is equal to r and eqs. 
( [42| ) through (5]) show that the spectral densities Sv, 
Si, and Svi satisfy the general relation (M), and since 
9|/W(6)| 2 /5e = -2Im{[/(+)(e)]*/(-)(e)}, the corre- 
lator Svi is also related to the response coefficient A by 
the expression identical to eq. (||). Moreover, since 

[/(+)( )]*j(-)(e) = -/|) + i / 1 / 2 sine, 







we see that the real part of the correlator Svi, which 
increases backaction dephasing produced by the SQUID, 
is indeed associated with the SQUID asymmetry. When 
Ii = h, the real part vanishes and the SQUID as detec- 
tor reaches quantum-limited optimum for measurement 
of the quantum coherent oscillations. For such a symmet- 
ric SQUID the signal-to-noise ratio of the oscillation mea- 
surement is given by eq. ( p9[ ) (with the intensity of back- 
ground output noise given by SV). If the temperature T 
of the symmetric SQUID is negligible, eqs. (42) through 
(H), and @ show that in this regime the SQUID is 
the quantum-limited detector with the energy sensitivity 
e = h/2, and the signal-to-noise ratio for measurement of 
the quantum flux oscillations is S max /Sv = 4. When T 
becomes non-vanishing, both the output and backaction 
noise increase, r' > r, and eq. (|2^) describes the gradual 
suppression of the signal-to-noise ratio with increasing 
temperature. 



B. Charge oscillations measured with a Cooper-pair 
electrometer 



Coherent oscillations of charge take place in Josephson 
junctions which are sufficiently small for the charging en- 
ergy Ec of an individual Cooper pair, E c = (2e) 2 /2C 
to be larger than temperature T and Josephson coupling 
energy Ej. The supercurrent flow through the junction 
is "discretized" in this regime into the transfer of individ- 
ual Cooper pairs by strong Coulomb repulsion. Quanti- 
tatively, if the charging energy is smaller that the super- 
conducting energy gap A so that the dissipative quasi- 
particle tunneling is suppressed, the junction dynamics 
is governed by the simple Hamiltonian: 

H = E c (n - qf - ^-(\n)(n + l\ + \n+ l)(n\) , (45) 

where n is the number of Coper-pairs charging the junc- 
tion, and, here and below, q is the charge (in units of 2e) 
injected into the junction from external circuit. Eigen- 
states of the Hamiltonian (|45|) form energy bands as func- 
tions of the injected charge q, which can be varied contin- 
uously. Variations of the injected charge q within these 
bands leads to the possibility of controlling the tunneling 
of individual Cooper pairsc3. 

The best way of injecting the charge q in a j|t 
tion is provided by the "Cooper-pair box" syster 
in which the junction is attached to external bias voltage 
V g through a capacitor. If q is fixed at half of a Cooper- 
pair charge, q ~ 1/2, and the tunneling amplitude Ej/2 
is much less than Ec, the two states of the Hamiltonian 
((45|): n — and n = 1 are nearly-degenerate and sepa- 
rated by the large energy gaps from all other states. In 
this regime, the junction dynamics is equivalent to that 
of a regular quantum two-state system with the two ba- 
sis state that correspond to a Coooper pair being on the 
left or on the right electrode of the junction. Coherent 



superposition of charge states in such a two-state system 
is observed indirectly by measuring either the width of 
the transition region between the two charge statesca or 
the energy gap between the eigenstatea!3 as functions of 
the induced charge q. Quantum coherent oscillations be- 
tween the two charge states were also observed directly 
in the time-dependent measuremcntlI3. 

C gl , I 



2e 



'0 



EjC, 
E 2 C 2 



R V, 



0/ 



FIG. 4. Schematic diagram of a continuous measurement of 
the quantum oscillations of a Cooper-pair in a small Joseph- 
son junction biased with charge V g C g ~ e through a small ca- 
pacitance C g . Capacitance Co couples the oscillating Cooper 
pair weakly to an electrometer composed of two junctions with 
Josephson coupling energies E\p and capacitances Ci,2- The 
electrometer is biased with a voltage Vb, while the tunneling 
current I through it is the measurement output. 

The experiment^ was effectively based on the strong 
measurement of charge oscillations, when each measure- 
ment suppresses the oscillations, and they are observed 
as oscillations of probability in an ensemble of measure- 
ments. Continuous weak measurement of the quantum- 
coherent charge oscillations similar to the measurement 
of the flux oscillations with a dc SQUID discussed above 
would provide a less intrusive way of studying these os- 
cillations. One of the detectors appropriate-ips such 
a measurement is a Cooper-pair electrometer! 3 HP !: two 
small Josephson junctions with Josephson coupling ener- 
gies E\ t 2 and capacitances C1.2 connected in series and 
shunted with a resistor R (Fig. 4). As in the case of dc 
SQUID, we consider dynamics of the Cooper pair trans- 
fer through the electrometer in the regime of incoherent 
tunneling. The main difference with the SQUID case 
is that now the amplitude of Cooper pair tunneling is 
modulated through the modulation of energy of the in- 
termediate state in the process of the two-step transfer of 
Cooper pairs in the two junctions. At small bias voltages 
Vb <C Ec/e, where from now on Ec — 2e 2 /(Ci + C2) is 
the charging energy of the central electrode of the elec- 
trometer, the intermediate state is virtual, and the aver- 
age current (/) through the electrometer is determined by 
the same eq. (^8|) with (O) replaced with the ampli- 
tude I(q) of Cooper-pair transfer through both junctions 
(defined below). The charge q injected into the central 
electrode controls the energy of the intermediate states 
in the process of the Cooper-pair transfer and modulates 
the tunneling amplitude. At small Josephson coupling 
energies E\j, -C Ec, and away from the resonance points 
q = ±1/2, Cooper-pair tunneling can be treated as per- 
turbation. Then the instantaneous value of the current 
I for a fixed value of the Josephson phase ip across the 
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electrometer (see, e.g., description of the two-junction 
system irH) is: 



I = I(q) sm(2eV t + ip(t)) , 



(46) 



Er 



l-2q l + 2q 



The two terms in the second equation in ( f46| ) correspond 
to the two intermediate states with different charges n = 
±1 on the central electrode of the electrometer in the 
Cooper-pair transfer process. 

Averaging eq. (^6|) over the equilibrium quantum fluc- 
tuations of ip we get expression for the average value of 
the current I that is equivalent to eq. (|38|). Since the 
tunneling current / through the electrometer is the mea- 
surement output, this expression determines the response 
coefficient A of the electrometer: 



X=(d(I)/dq)=7r(d[I(q)] 2 /dq)T/e : 



(47) 



where r is given by eqs. (|3^) and (f40|). Similarly, the 
output noise of the current I can be obtained as 



(48) 



The backaction noise of the Cooper-pair electrometer 
is created by fluctuations of the charge on its central elec- 
trode in the process of the Cooper-pair tunneling. This 
fluctuations lead to fluctuations of electric potential of 
this electrode. In the same regime as for eq. (|46|), the 
magnitude of this fluctuations is determined by the mag- 
nitude of the instantaneous value of the potential at a 
fixed Josephson phase difference <p across the electrome- 
ter: 



V = U(q) cos(2eV t + <p(t)) , 



(49) 



where 



U(q) 



2eE c y (l~2q) 2 (1 + 2q) 



As before, averaging over the fluctuations of (p(t) we get 
the spectral density of the backaction noise and its cor- 
relation with the output noise: 



5', 



Mq)] 



V. 



Svi = -iU{q)I(q)T' . 



(50) 



Equations (fl7|), (|||), and ([so]) show that at vanishing 
temperature, when r' = r, the noise characteristics of 
the Cooper-pair electrometer satisfy the general relations 
(fjj) and (||) of a quantum-limited detector. Moreover, the 
electrometer is "symmetric" detector in a sense that the 
input-output correlator Syi is purely imaginary. This 
means that both its energy sensitivity e and the signal-to- 
noise ratio ( p9| ) for measurement of the two-state system 
reach fundamental limits. 

In summary, the examples of specific detectors consid- 
ered in this Section show explicitly that many standard 



detectors should be capable of reaching the fundamental 
limits of sensitivity for measurements of electric charge 
and magnetic flux. In this regime, they are characterized 
by the fundamental signal-to-noise ration of 4 for con- 
tinuous weak measurement of the macroscopic quantum 
coherent oscillations. The limitation on the signal-to- 
noise ratio of such a measurement has the same origin 
as the quantum limitation on the operation of an ideal 
linear phase-insensitive amplifier that adds a minimum 
of half-a-quantum of noise to the amplified signal. 
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Likharev, J.E. Lukens, Yu.V. Nazarov, R.J. Schoelkopf, 
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Appendix 

In this Appendix, we derive eq. ([32]) that relates the 
matrix elements of the perturbation of the scattering po- 
tential to the transmission properties of a point contact. 
Relation ([32]) can be established considering the station- 
ary states of an electron confined to move on the interval 
x 6 [— L/2, L/2] with the main scattering potential lo- 
cated at the center of the interval, x ~ 0. The scattering 
matrix S for the symmetric scattering potential can be 
written as 



S = 



i sin v . 



cos v 
i sin v 



where O is the phase of the transmission amplitude, and 
v parametrizes transmission probability D: 



D 



cos 2 v . 



(51) 



Diagonalizing the scattering matrix S, we find the two 
phase shifts (fi^ associated with it: ipi = + v, <pa = 
O + 7r — v. The variations Sipj of the two phase shifts 
due to perturbation of the scattering potential lead to 
changes 5ej in energies of the two stationary states, 



5sj — VF^tfij/L 



(52) 



For symmetric main scattering potential, the two sta- 
tionary states Xj ( x ) are given by the even and odd com- 
binations of the scattering states ipj (x) . In the basis of 
the states Xji the perturbation matrix U introduced in 
eq. (|0]) is: 



U 



(U n + U 22 )/2 + RcU 12 , (ETu - U 22 )/2 - ilmU 12 
(U n - U 22 )/2 + ilmU 12 , (U u + U 22 )/2 - RcU 12 

(53) 



The diagonal elements of this matrix give the first-order 
corrections to the energies of the stationary states. Com- 
paring expressions for the energy corrections given by eq. 
(53) to expressions for the phase shifts combined with eq. 
(52), we see that 
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Ret/i 



This equation, together with the relation (|5l| ) between v 
and transmission probability D, gives eq. (p2f) of the main 
text. Since the matrix ( j53|) should be symmetric for the 
perturbation of the scattering potential that is symmetric 
with respect to the main part of the potential, we also see 
that Im[/i2 = in the symmetric case. This means that 
the nonvanishing imaginary part of U12 can be viewed 
as a measure of the asymmetry of coupling between the 
point contact and a source of the perturbation. 
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